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Abstract 

We study the decay of D'^ and Df mesons into charged five-body final states, and 
report the discovery of the decay mode K~^K~'k~ , as well as measure¬ 

ments of the branching ratios of the decay modes D'^ K~K A Dt - 
K^K TT'^TT'^TT , Df (pTr'^TT'^TT and D'^/Df 7r'*'7r'*'7r'*“7r tt . An analysis of 
the resonant substructure for ^ A“7r'''7r'''7r'''7r“ and Df -i- K~^K~n~^n~^n~ is 
also included, with evidence suggesting that both decays proceed primarily through 
an ai vector resonance. 

PACS numbers; 13.25.Ft, 14.40Lb 


The hadronic hve-body decays of charmed mesons have been studied in recent 
years [1-6], but limited statistics have prevented precise measurements of their 
resonant substructure. Theoretical predictions are limited mainly to two-body 
decay modes, and little is known about how hve-body hnal states are produced. 
Theoretical discussion suggests a “vector-dominance model,” in which heavy 
havor mesons decay into a two-body intermediate state by emitting a W, which 
immediately hadronizes into a charged vector, axial vector, or pseudoscalar 
meson [7]. The charged meson then decays strongly to produce a many-body 
hnal state. Conhrmation of this model could provide a mechanism for the 
production of hve-body hnal states. 

The FOCUS Collaboration [8-10] has studied two hve-body decay modes. 
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K~ Ti^Ti^Ti~ and K'^K~n~^'K^'K~. We find evidence that in 

both modes the resonant substructure is dominated by a two-body vector 
resonance involving the ai(1260)’''. We also present inclusive branching ratio 
measurements of four charged five-body hadronic decays, including the first 
evidence of the decay mode ^ K~'K^7i^n~. 

Five-body D+ and D+ decays are reconstructed using a candidate driven ver¬ 
tex algorithm [ 8 ]. A decay vertex is formed from the hve reconstructed tracks. 
The momentum vector of the parent D meson is then used as a seed to inter¬ 
sect other tracks in order to hnd the production vertex. Events are selected 
based on a number of criteria. The confidence level of the decay vertex is 
required to be greater than 1%. The confidence level that a track from the 
decay vertex intersects the production vertex is required to be less than 1 %. 
The likelihood for each particle to be a proton, kaon, pion, or electron based 
on Cerenkov particle identihcation is used to make additional requirements [9]. 
For each kaon candidate we require the negative log-likelihood kaon hypoth¬ 
esis, Wk = —2 ln(kaon likelihood), to be favored over the corresponding pion 
hypothesis ITC by 14A — Wk > 3. In addition, for each pion candidate we 
require the pion hypothesis to be favored over any alternative hypothesis. We 
also require the signihcance of separation of the production and decay vertices 
to be at least 10. In order to reduce background due to secondary interactions 
of particles from the production vertex, we require the D reconstructed mo¬ 
mentum to be greater than 25 GeV/c and the secondary vertex to be outside 
of the target. Finally, we remove events that are consistent with various D* 
decays. 

We turn now to additional analysis cuts made in individual modes, beginning 
with D~^ K~7r~^7r~^7i~^7r~ . Because this mode is the most abundant we apply 

only the standard cuts used in all modes. Figure la shows the KAti invariant 
mass plot. The distribution is fitted with a Gaussian for the D~^ signal and a 
2 °d degree polynomial for the background. A binned maximum likelihood £t 
gives 2923 ± 78 events. 

The D~^/D^ 7r~^7r~^7r~^7r~7r~ modes are more difficult to detect, due to the 
large combinatorial background. To reduce this background we increase the 
separation of the secondary vertex from the target to two standard deviations. 
We further impose a series of selection cuts to remove misidentihed charm 
decays. We remove the decays ^ ri'7r~^, 77 '—> 7 r+ 7 r+ 7 r“ 7 r“ 7 r° by requiring 

the four pion reconstructed mass to be larger than the 77 ' — 7 r° mass difference, 
that is, M 4 ^ > 0.825 GeV/c^. Figure lb shows the hve-pion invariant mass plot 
for events that satisfy these cuts. The distribution is fitted with a Gaussian for 
the D~^ signal (835 ±49 events), another Gaussian for the Df signal (671 ±47 
events) and a 1 ®* degree polynomial for the background. 

For the K~ t:~ mode the requirement of two kaons in the final 
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state greatly reduces background, allowing us to apply only the standard cuts 
used in all modes. Figure Ic shows the K~ti~ invariant mass plot for 
events satisfying these cuts. We fit to a Gaussian (240 ± 30 events) and 2°'^ 
degree polynomial. 

For the hnal state we have also studied the subresonant decay 

'K~, by additionally requiring the K^K~ invariant mass combi¬ 
nation to be consistent with the 0 mass. The (j)'K^'K^'K~ invariant mass plot 
is shown in Figure le. We £t to a Gaussian (136 ± 14 events) and 2°'^ degree 
polynomial. 

The decay ^ K~^K~n~^n~^n~ is Gabibbo suppressed. We require a signif¬ 
icance of vertex separation of 20, and tighten particle identification cuts on 
both kaons to — Wk > 4, but remove all requirements on the pions. We 
also require the reconstructed momentum to be greater than 50 GeV/c. 
Figure Id shows the resulting K~'k~ invariant mass plot. This is the 
first observation of this mode. We fit with a Gaussian for the signal (38 ±8 
events), another Gaussian for the Dg events and a 2^^^ degree polynomial for 
the background. 

We measnre the branching fraction of the K~t[~ mode rela¬ 
tive to K~7r~^7r~^, then measnre the branching fractions of the other 

D~^ modes relative to the D'^ K~n~^7r'^7r~^7r~ to rednce systematic effects 

dne to differences in the nnmber of decay prodncts. All decay modes are 
measured relative to Df K^K~'k^ . The normalizing decay modes are snb- 
jected to the same vertex cnts and analogous Cerenkov identification cuts as 
the mode in qnestion to minimize systematic errors. The detector and analysis 
efficiency is calculated using a Monte Carlo simulation. For modes included in 
onr resonant substructure analysis the Monte Carlo contains the incoherent 
mixtnre of snbresonant decays determined by our analysis. For modes not in- 
clnded in onr resonant snbstructnre analysis, the Monte Carlo is composed of 
five-body phase space. We test for dependency on cnt selection by individn- 
ally varying each cnt. The results, compared with existing measnrements, are 
shown in Table 1. 

We stndied systematic effects dne to nncertainties in Monte Carlo simulation, 
fitting procednre, resonant substrncture, Monte Carlo statistics and absolnte 
tracking efficiency. To determine the systematic error we follow a procednre 
based on the S-factor method used by the Particle Data Group [11]. For each 
mode we split the data sample into fonr independent snbsamples based on D 
momentnm and period of time in which the data was collected. We then define 
the split sample variance as the difference between the scaled variance and the 
statistical variance if the former exceeds the latter. We also evalnate systematic 
effects associated with Monte Carlo simnlation of multi-body decays. The 
branching ratios are evalnated with mnltiple conditions on the isolation of the 
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Fig. 1. (a) Kin invariant mass distribution, (b) Svr invariant mass distribution, (c) 
KK^tt invariant mass distribution for optimized cuts, (d) KKStt invariant mass 
distribution for D~^ optimized cuts, (e) (pSir invariant mass distribution.The fits are 
described in the text and the numbers quoted are the yields. 

production vertex, with the variance used as the systematic error. In addition 
we evaluate the systematic effects based on different htting procedures. The 
branching ratios are evaluated under various £t conditions, and the variance 
is used as the systematic error, as all fit variants are a priori equally likely. We 
also evaluate systematic effects due to uncertainty in resonant substructure by 
calculating the branching ratios using various mixtures of subresonant states in 
the Monte Carlo. The variance in the branching ratios for different subresonant 
mixtures is used as the systematic error, treating each subresonant mixture as 
a priori equally likely. We also evaluate the systematic effect from Monte Carlo 
statistics, adding in quadrature the uncertainty in the calculated efficiencies 
from both the signal and normalizing mode. Finally we evaluate systematic 
effects from uncertainty in absolute tracking efficiency of multi-body decays 
using studies of ^ K~TT~^'n'~^TT~ and K~'k^ decays. The systematic 

effects are then all added together in quadrature to obtain the final systematic 
error. 

In addition to reporting inclusive branching ratio measurements, we have stud- 
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Table 1 

Branching ratios for five-body modes and comparison to the previous measurements 
by E687. All branching ratios are inclusive of subresonant modes. 


Decay Mode 

FOCUS 

E687[6] 

r ( D ~ TT TT ■*" TT ■*" TT ~ ) 

0.058±0.002±0.006 

0.290±0.017±0.011 

0.145±0.011±0.010 

0.150±0.019±0.025 

0.249±0.024±0.021 

0.040±0.009±0.019 

0.077±0.008±0.010 

0.299±0.061±0.026 

0.158±0.042±0.031 

0.188±0.036±0.040 

0.28±0.06±0.01 

r(D+^7^-7r+7r+) 

r { Z? TT TT TT TT “ TT ~ ) 

77+77 + 77*^77”) 
r { Z? > 77 + 77 + 77 + 77 ~ 77 ~ ) 

V{Dt^K-K+-K+) 

T{Dt^K+K--K+-K+-K-) 

r(D+^i^-x+77+) 

r ( Z? (/) 77 + 77 + 77 ~ ) 

r{D+^4,-K+) 

T{D+^K+K--w+tt+tt-) 

rfZ)+— >ZC~77 + 77 + 77 + 77”) 



ied the resonance substructure in two decays: ^ K~TT~^n^n^n~ and 

K~7r'^TT~^7r~. We use an incoherent binned fit method, also used in reference 
6 , a simplihed approach which assumes the hnal state is an incoherent super¬ 
position of subresonant decay modes containing vector resonances. A coherent 
analysis for decays into hve-body hnal states has not yet been attempted, and 
would be very difficult given the statistics of this experiment. For the —>■ 
K~7r~^7r^7r^7r~ mode we consider the lowest mass {K~n^) and {n^n~) reso¬ 
nances, as well as a nonresonant channel; K*^7r~7r~^n~^, K~p^7v~^7i~^, K*^p^7r~^, 
and (iF“7r’''7r’''7r’''7r“)NR- All states not explicitly considered are assumed to be 
included in the nonresonant channel. 

We determine the acceptance corrected yield into each subresonant mode us¬ 
ing a weighting technique whereby each event is weighted by its kinematic 
values in three submasses: ( 7 r’'' 7 r“), and ( 7 r’'' 7 r’''). No resonance in the 

submass exists, but we include it in order to compute a meaningful 
estimate of the ht. The weights are obtained using separate Monte Carlo sim¬ 
ulations for the four decay modes, with the Particle Data Group values [11] for 
the mass and width of each resonance. Eight kinematic bins are constructed 
depending on whether each of the three submasses falls within the expected 
resonance (In the case of 7 r+ 7 r+, the bin is split into high and low mass re¬ 
gions). For each Monte Carlo simulation the bin population in the eight bins is 
determined using a sideband subtracted cut on the peak, allowing a linear 
transformation matrix to be calculated. The weights are then determined from 
the transformation matrix by a minimization procedure. Each data event 
which satishes our selection cuts is then weighted according to its kinematic 
values in the submass bins. Once the weighted distributions for each of the 
four modes are generated, we determine the acceptance corrected yield by £t- 
ting the distributions with a Gaussian signal and a linear background. Using 
incoherent Monte Carlo mixtures of the four subresonant modes we verified 
that our procedure was able to correctly recover the generated mixtures of the 
four modes. 
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Fig. 2. K TT~^Tr~^Tr~^TT weighted invariant mass for (a) {K 7r~^Tr~^7r ~^tt )nr, (b) 
K*^7r~7:~^7:~^, (c) K~, (d) K*^p^7:~^, (e) Inclusive sum of all four modes. 

The results for K~TT~^7r'^7r~^7r~ are summarized and compared to the E687 re¬ 
sults in Table 2. Taking into account the correlation among the subresonant 
fractions, the calculated fh® hypothesis that the results are consistent 

with E687 is 6.5 (4 degrees of freedom). The four weighted histograms with 
hts are shown in Fig. 2. Figure 2e is the weighted distribution for the sum of 
all subresonant modes. The goodness of £t is evaluated by calculating a 
for the hypothesis of consistency between the model predictions and observed 
data yields in each of the 8 submass bins. The calculated x^ is 7.4 (4 degrees 
of freedom), with most of the contribution resulting from a poor Monte 
Carlo simulation of the spectrum for the mode. We assessed 

systematic errors by individually varying the width of the submass bins cor¬ 
responding to the p and resonances by 20%. The systematic error is then 
estimated as the variance of the two measurements with varied widths, along 
with the original measurement. Since our methods of calculating subresonant 
fractions and inclusive branching ratios are distinct, statistical and systematic 
errors are added in quadrature when normalizing our subresonant fractions to 
other modes. 

We follow a similar procedure for the , treating the 
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Table 2 

Fractions relative to the inclusive mode and comparison to previous measurements 
for the resonance substructure of the —> iF“7r^7r^7r^7r“ decay mode. These 

values are not corrected for unseen decay modes. 


Subresonant Mode 

Fraction of Kdvr 

E687 Fraction [6] 

(iF“ vr"'' TT"''TT^ 7r“ ) nr 

0.07T0.05T0.01 

< 0.26(90%C.L.) 

K*0j^-^+^+ 

0.21T0.04T0.06 

0.42T0.14 

K~ p^t:~^7:~^ 

0.30T0.04T0.01 

0.44T0.14 

K*^p^7r~^ 

0.40T0.03T0.06 

0.20T0.09 


final state as an incoherent superposition of the {K^K~) and ( 7 r’'' 7 r“) reso¬ 
nances, as well as a nonresonant channel: K~, (ppTr~^ and 

(ii'+K“7r+7r+7r“)NR- Each event is weighted by its value in each of three sub¬ 
masses: and and the weighted distributions are 

again htted with a Gaussian signal and a linear background. The results are 
summarized in Table 3 and are presented in Fig 3. The goodness of fit is 
evaluated by calculating a for the hypothesis of consistency between the 
model predictions and observed data yields in each of the eight submass bins. 
The calculated is 10.2 (4 degrees of freedom), with most of the x^ contri¬ 
bution resulting from a poor Monte Carlo simulation of the spectrum 

in the nonresonant channel. We assess systematic errors by calculating the 
variance of our results with 20 % variations in the width of the submass bins 
corresponding to the p and 0 resonances. 

In both resonant substructure analyses the dominant mode is of the form 
vector-vector-pseudoscalar: K*'^p^Ti^ and (j)p^Ti^ in the case of 
and K~'K~, respectively. Given the phase space constraints for both 
of these decays, such a result is unexpected. However, theoretical discussion 
of a vector-dominance model for heavy flavor decays [7] suggests that charm 
decays are dominated by quasi-two-body decays in which the immediately 
hadronizes into a charged pseudoscalar, vector, or axial vector meson. Thus 
branching ratios of the form D ai(1260)’''W are of comparable value to 

Table 3 

Fractions relative to the inclusive mode for the resonance substructure of the D'^ 
K ~^decay mode. These values are not corrected for unseen decay modes. 


Subresonant Mode 

Fraction of 2K37r 

( 77 77 “ TT ■*“ TT “*■ TT “ ) N R 

0.10T0.06T0.05 

07r“7r'’'7r'’' 

0.21T0.05T0.06 

77+i7-p°7r+ 

<0.03 (90% C.L.) 

(l)p^ 7r~^ 

0.75T0.06T0.04 
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Fig. 3. K^K weighted invariant mass for (a) {K^K )nr, (b) 

(/)7r“7r+7r+, (c) K~p^ir^, (d) (j)p^7T^, (e) Inclusive sum of all four modes. 

those observed for n^X, when adjusted for phase space. Such theoretical 
discussion raises the possibility that the resonant substructure for both modes 
is dominated by a quasi-two-body decay involving the ai; K*^ai and (paf for 
K“7r+7r+7r+7r“ and i^+K“7r+7r+7r“, respectively, where af p^7r~^. Although 
the central value of the ai mass lies outside of phase space for both decays, 
these decay modes are allowed due to the large width of the ai. However, 
the large width of the ai and its position in phase space make the resonance 
difficult to detect directly. 

To verify the subresonant decays are proceeding through ai we generate Monte 
Carlo simulations of —>■ K*^af and Df (pa-t^ assuming the Oi has a width 
of 400 MeV/c^ and decays entirely as an S-wave, and use our subresonant 
analysis procedure explained above. In both cases the yield fractions in each 
of the subresonant modes from the oi Monte Carlo are similar to the reported 
fractions from the data, with particular agreement in the case of Df —> (paf. 
Such agreement suggests both channels may be dominated by a two-body 
intermediate state involving the ai. 

Accepting the hypothesis that hve-body modes are dominated by quasi-two- 
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body decays, we calculate branching ratios for the decays K*^ai and 

(pOi using the ratios of the observed fractions of K *^and (ppn from 
data (40% and 75%) to the observed fractions from Monte Carlo simulations 
of K*^ai and D+ ^ (j)ai (70% and 78%). Assuming the af decays to 

p^7r~^ 50% of the time and using the Particle Data Group 0 and K*^ branching 
fractions [11], the D~^ ^ K*^ai and Df (paf branching fractions, including 
unseen decays, are shown in Table 4. We assess systematic errors by increasing 
the width of the ai resonance in our generated Monte Carlo to 600 MeV/c^, 
taking the systematic error as the variance of our measurements with the two 
widths. 

Table 4 

Inclusive branching ratios for states. These values are corrected for unseen decay 
modes. _ 


Decay Mode 

Fraction 

T{D+^K*°a+) 

r(D+^iC-7r+7r+) 

r(D+^flia+) 

0.099±0.008±0.018 

0.559±0.078±0.044 

TiDt^K+K--K+) 


In conclusion we have measured the relative branching ratios of hve-body and 
three-body charged hadronic decays of D~^ and and have presented the hrst 
evidence of the decay mode D~^ K~'k~ . We have also performed an 

analysis of the resonant substructure in the decays and 

K~ 'K~. Our analysis provides some evidence that both decays 
proceed through a quasi-two-body decay involving the ai(1260)’'' particle. 
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